約 46,524 件
https://w.atwiki.jp/ce00582/pages/1975.html
public class tax17{ public static void main(String[] args) { double th1[]=new double[11]; double th2[]=new double[11]; double y1[][]=new double[11][11]; double y2[][]=new double[11][11]; double c[][]=new double[11][11]; double tk,tl,tr; double c1,w1,w2,l1,l2; double we; double maxwel,maxtk,maxtl,maxtr; int n1,n2,s1,s2; for (s1=1;s1 11;s1++){ th1[s1]=0.2*s1; } for (s2=1;s2 11;s2++){ th2[s2]=0.1*s2; } maxwel=-999; maxtk=0; maxtl=0; maxtr=0; for (n1=5;n1 40;n1++){ for (n2=5;n2 40;n2++){ tk=0.01*n1; tl=0.01*n2; tr=trs(tk,tl,th1,th2); we=wel(tk,tl,tr,th1,th2);; if (we maxwel)maxtk=tk; if (we maxwel)maxtl=tl; if (we maxwel)maxtr=tr; if (we maxwel)maxwel=we; } } tk=maxtk; tl=maxtl; tr=maxtr; for (s1=1;s1 11;s1++){ for (s2=1;s2 11;s2++){ l1=mlx(s1,s2,tk,tl,tr,th1,th2); l2=flx(s1,s2,tk,tl,tr,th1,th2); w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; y1[s1][s2]=th1[s1]*l1; y2[s1][s2]=th2[s2]*l2; c[s1][s2]=w1*l1+w2*l2+tr; } } System.out.println(tk); System.out.println(tl); } static double bud(double tk,double tl,double tr,double th1[],double th2[]){ double b1,w1,w2,l1,l2,c1; int s1,s2,i; b1=0; for (s1=1;s1 11;s1++){ for (s2=1;s2 11;s2++){ w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; l1=1-(w1+w2+tr)/(3*w1); l2=1-(w1+w2+tr)/(3*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=0; if (i==2)l2=(w2-tr)/(2*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=(w1-tr)/(2*w1); if (i==2)l2=0; c1=w1*l1+w2*l2+tr; b1=b1+th1[s1]*l1+th2[s2]*l2-c1; } } return b1; } static double wel(double tk,double tl,double tr,double th1[],double th2[]){ double we,w1,w2,l1,l2,c1; int s1,s2,i; we=0; for (s1=1;s1 11;s1++){ for (s2=1;s2 11;s2++){ w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; l1=1-(w1+w2+tr)/(3*w1); l2=1-(w1+w2+tr)/(3*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=0; if (i==2)l2=(w2-tr)/(2*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=(w1-tr)/(2*w1); if (i==2)l2=0; c1=w1*l1+w2*l2+tr; we=we+Math.log(c1)+Math.log(1-l1)+Math.log(1-l2); } } return we; } static double trs(double tk,double tl,double[] th1,double[] th2){ double tr1,tr2,tr3,b1,b2; int t1; tr1=0.1; tr2=0.05; b1=bud(tk,tl,tr1,th1,th2); t1=0; while(t1 100){ b2=bud(tk,tl,tr2,th1,th2); tr3=tr2-b2*(tr2-tr1)/(b2-b1); b1=b2; tr1=tr2; tr2=tr3; if (b2*b2 0.001)t1=1000; t1=t1+1; } return tr2; } static double ux(double c1,double l1,double l2){ int pp; double u1,cp1,lp1,lp2; pp=0; cp1=c1; lp1=l1; lp2=l2; if (lp1 0)pp=100; if (lp1 0.99)pp=100; if (lp2 0)pp=100; if (lp2 0.99)pp=100; if (pp 50)lp1=0.5; if (pp 50)lp2=0.5; u1=Math.log(cp1)+Math.log(1-lp1)+Math.log(1-lp2); if (pp 50)u1=-999; return u1; } static double mlx(int s1,int s2,double tk,double tl,double tr,double[] th1,double th2[]){ int i; double w1,w2,l1,l2; w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; l1=1-(w1+w2+tr)/(3*w1); l2=1-(w1+w2+tr)/(3*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=0; if (i==2)l2=(w2-tr)/(2*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=(w1-tr)/(2*w1); if (i==2)l2=0; return l1; } static double flx(int s1,int s2,double tk,double tl,double tr,double[] th1,double th2[]){ int i; double w1,w2,l1,l2; w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; l1=1-(w1+w2+tr)/(3*w1); l2=1-(w1+w2+tr)/(3*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=0; if (i==2)l2=(w2-tr)/(2*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=(w1-tr)/(2*w1); if (i==2)l2=0; return l2; } }
https://w.atwiki.jp/yaruobetabeta/pages/29.html
l二二"'i [][]`ヽ / / | !-──`ヽ; ',...,,_ ,.. '"´ヽ、__ | | ' ; `ヽ.ノ / ゴ , '"'; ' , '; ! `',. ' ; -'`ヽ; `ヽ, l二二..,」 \ く ,' ゝ、_ `ヽ. i/ / /`' 、,i i r─ ┐[][]. ', `ヽ; ', /-、 / / ,' ;'、 二l ! rー'"7´`''ー--、!__,,..イゝ、 >i ハイ<., / ; イ ヽ. └─┘,ノヽ、_/ `ヽ.,__,. '"´ ̄`ヽ、`ヽ `ヽ< / i /! '; !_ `ヽr'"´`>r-‐'´ヽ. `7、\ `ヽ、.,____,,/ヽ、__/ | ヽ、 `ヽ/ .,. '´ .ハ 丶rァ'"`ヽ、 / `7´ ̄`ヽ; r''"`ヽ、 ノ / ヽ,.' / / ! ; L..__ノ、 `Y`ヽ、 `7´ / / ,' /`ヽ!, ,ハ i `ヽ. /! `'r-'、イ i .,' ,ァ'、 !`メ i ,i ! i Y`ヽ! ,' ト、 ', i イ i`>'、、, | /! ハ ,ハ ハ ', ヽ、__,.ヘ. / i ! ノ ! !.'、 '、 ,り`'レ' ! /,、!,_/-‐'ァ' / !. ', Y .| ∠..__ノ| /!. ` `ニ´ レ' !`"ァーr'、,' / ! i⌒ヽ! レ' ! "''' '、´ り ノ! ,.' ; i | ノ ,/ 人 , `ニ´' /イ /i ,' |イ / ,' / .ノ`ヽ., ,.,.,. ,. ' .レ' | / ', ヽ、 ヽ ,' i/ ( ノヽ、 ` ー-‐ / .l二二"'i [][] i ' , ァ、ヽ <O ノ ノ!_,,..ソ-─'-r- 、,ヽ.,,_______,,.. ,,.イ,' ;. | | ノ ノ、 ノ i `'' ,.( 〈´ ', i `''<´ // ハ. l二二..,」 r─ ┐[][] `ヽ. / ;' '; ノ '"' , ',. ! / ヽ.イ .i. / i i ト、 二l ! ' , i i !O> / i / ' , ' , ノ'´) / ヽ. ', ! / ,' ,' └─┘ ゴ `ヽ、 i '; ', ,' / ', rイ ヽ ヽ,.' ,.'ァ´`i-、 / ヽ.Y / / ,.イ `ヽ; `ヽ i i / / <O>ノ ト、 iヽ/ i/ / i、' ,. ヽi,.イ / ! ヽ; Y! ',/ <O> 八雲紫(東方Project) 初出はEX1スレ目2379 2年ぶりに再開したおっぱい幼馴染 ちなみにちっぱい幼馴染だった場合は佐々木であった 安価に嫌われていたが、EX1スレ5503にてようやくデートの約束に至る、すすめゆかりんおっぱいおっぱい 趣味は2c(ry パソコン、動いたりするのは好きじゃない EXVSでアッガイで50人抜き達成 Q、友達は、いますか? A,やる夫がいるじゃない バストサイズはCカップくらい改め、高校生にしちゃ大きいくらいである
https://w.atwiki.jp/javadsge/pages/8774.html
import numpy as np import matplotlib.pyplot as plt x = np.random.randn(30) y = np.sin(x) plt.plot(x, y, "go") plt.show()
https://w.atwiki.jp/ce00582/pages/2018.html
public class tax111{ public static void main(String[] args) { double th1[]=new double[11]; double th2[]=new double[11]; double y1[][]=new double[11][11]; double y2[][]=new double[11][11]; double c[][]=new double[11][11]; double up[][][][][]=new double[11][11][3][3][3]; double wp[][][][][]=new double[11][10][3][3][3]; int opc[][][]=new int[11][11][101]; int op1[][][]=new int[11][11][101]; int op2[][][]=new int[11][11][101]; double tk,tl,tr; double c1,w1,w2,l1,l2,u1,u2; double we; double maxwel,maxtk,maxtl,maxtr; double h; int n,n1,n2,n3,s1,s2,nx1,nx2,nx3; int t1,t2,t3; int j,ep; for (s1=1;s1 11;s1++){ th1[s1]=0.2*s1; } for (s2=1;s2 11;s2++){ th2[s2]=0.1*s2; } maxwel=-999; maxtk=0; maxtl=0; maxtr=0; for (n1=5;n1 40;n1++){ for (n2=5;n2 40;n2++){ tk=0.01*n1; tl=0.01*n2; tr=trs(tk,tl,th1,th2); we=wel(tk,tl,tr,th1,th2);; if (we maxwel)maxtk=tk; if (we maxwel)maxtl=tl; if (we maxwel)maxtr=tr; if (we maxwel)maxwel=we; } } tk=maxtk; tl=maxtl; tr=maxtr; for (s1=1;s1 11;s1++){ for (s2=1;s2 11;s2++){ l1=mlx(s1,s2,tk,tl,tr,th1,th2); l2=flx(s1,s2,tk,tl,tr,th1,th2); w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; y1[s1][s2]=th1[s1]*l1; y2[s1][s2]=th2[s2]*l2; c[s1][s2]=w1*l1+w2*l2+tr; } } h=0.001; for (s1=1;s1 11;s1++){ for (s2=1;s2 11;s2++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (n3=-1;n3 2;n3++){ c1=c[s1][s2]+n3*h; l1=(y1[s1][s2]+n1*h)/th1[s1]; l2=(y2[s1][s2]+n2*h)/th2[s2]; up[s1][s2][n1+1][n2+1][n3+1]=ux(c1,l1,l2); } } } } } for (s1=1;s1 11;s1++){ for (s2=1;s2 10;s2++){ for (n1=-1;n1 2;n1++){ for (n2=-1;n2 2;n2++){ for (n3=-1;n3 2;n3++){ c1=c[s1][s2]+n3*h; l1=(y1[s1][s2]+n1*h)/th1[s1]; l2=(y2[s1][s2]+n2*h)/th2[s2+1]; wp[s1][s2][n1+1][n2+1][n3+1]=ux(c1,l1,l2); } } } } } for (s1=1;s1 11;s1++){ for (n=1;n 101;n++){ t3=0; while(t3 100){ n1=0; n2=0; n3=0; t1=0; s2=1; while (t1 100){ n1=chance(); n2=chance(); n3=chance(); u1=up[s1][s2][n1+1][n2+1][n3+1]; if (u1 -100)t1=1000; } opc[s1][s2][n]=n3; op1[s1][s2][n]=n1; op2[s1][s2][n]=n2; for (s2=2;s2 11;s2++){ t2=0; n1=0; n2=0; n3=0; while (t2 100){ n1=chance(); n2=chance(); n3=chance(); u1=up[s1][s2][n1+1][n2+1][n3+1]; j=0; if (u1 -100)j=1000; nx1=op1[s1][s2-1][n]; nx2=op2[s1][s2-1][n]; nx3=opc[s1][s2-1][n]; w1=wp[s1][s2-1][nx1+1][nx2+1][nx3+1]; if (w1 u1)j=0;; if (j 100)t2=1000; } opc[s1][s2][n]=n3; op1[s1][s2][n]=n1; op2[s1][s2][n]=n2; } ep=0; for (s2=1;s2 11;s2++){ ep=ep+op1[s1][s2][n]+op2[s1][s2][n]-opc[s1][s2][n]; } j=0; if (ep 2)j=j+1; if (ep -2)j=j+1; if (j==2)t3=1000; } } } } static int chance(){ double x; int xp; x=Math.random(); xp=0; if (x 0.7)xp=1; if (x 0.7)xp=-1; return xp; } static double bud(double tk,double tl,double tr,double th1[],double th2[]){ double b1,w1,w2,l1,l2,c1; int s1,s2,i; b1=0; for (s1=1;s1 11;s1++){ for (s2=1;s2 11;s2++){ w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; l1=1-(w1+w2+tr)/(3*w1); l2=1-(w1+w2+tr)/(3*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=0; if (i==2)l2=(w2-tr)/(2*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=(w1-tr)/(2*w1); if (i==2)l2=0; c1=w1*l1+w2*l2+tr; b1=b1+th1[s1]*l1+th2[s2]*l2-c1; } } return b1; } static double wel(double tk,double tl,double tr,double th1[],double th2[]){ double we,w1,w2,l1,l2,c1; int s1,s2,i; we=0; for (s1=1;s1 11;s1++){ for (s2=1;s2 11;s2++){ w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; l1=1-(w1+w2+tr)/(3*w1); l2=1-(w1+w2+tr)/(3*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=0; if (i==2)l2=(w2-tr)/(2*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=(w1-tr)/(2*w1); if (i==2)l2=0; c1=w1*l1+w2*l2+tr; we=we+Math.log(c1)+Math.log(1-l1)+Math.log(1-l2); } } return we; } static double trs(double tk,double tl,double[] th1,double[] th2){ double tr1,tr2,tr3,b1,b2; int t1; tr1=0.1; tr2=0.05; b1=bud(tk,tl,tr1,th1,th2); t1=0; while(t1 100){ b2=bud(tk,tl,tr2,th1,th2); tr3=tr2-b2*(tr2-tr1)/(b2-b1); b1=b2; tr1=tr2; tr2=tr3; if (b2*b2 0.001)t1=1000; t1=t1+1; } return tr2; } static double ux(double c1,double l1,double l2){ int pp; double u1,cp1,lp1,lp2; pp=0; cp1=c1; lp1=l1; lp2=l2; if (lp1 0)pp=100; if (lp1 0.99)pp=100; if (lp2 0)pp=100; if (lp2 0.99)pp=100; if (pp 50)lp1=0.5; if (pp 50)lp2=0.5; u1=Math.log(cp1)+Math.log(1-lp1)+Math.log(1-lp2); if (pp 50)u1=-999; return u1; } static double mlx(int s1,int s2,double tk,double tl,double tr,double[] th1,double th2[]){ int i; double w1,w2,l1,l2; w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; l1=1-(w1+w2+tr)/(3*w1); l2=1-(w1+w2+tr)/(3*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=0; if (i==2)l2=(w2-tr)/(2*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=(w1-tr)/(2*w1); if (i==2)l2=0; return l1; } static double flx(int s1,int s2,double tk,double tl,double tr,double[] th1,double th2[]){ int i; double w1,w2,l1,l2; w1=(1-tk)*th1[s1]; w2=(1-tl)*th2[s2]; l1=1-(w1+w2+tr)/(3*w1); l2=1-(w1+w2+tr)/(3*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=0; if (i==2)l2=(w2-tr)/(2*w2); i=0; if (l1 0)i=i+1; if (l2 0)i=i+1; if (i==2)l1=(w1-tr)/(2*w1); if (i==2)l2=0; return l2; }
https://w.atwiki.jp/1548908-wdt1/pages/467.html
パック デュエル・ターミナル EX1(P)DT1:DUEL TERMINAL -シンクロ覚醒!!-(OCG):DUEL TERMINAL -ワームの侵攻!!-(OCG) BEFORE:STARSTRIKE BLAST(P)DT1 NEXT デュエル・ターミナル EX2(P)DT1 条件 滞在時:ビフレスト01F・04F 滞在時:フリズスキャルヴ・ニヴルヘイム02F 開放鍵:ニヴルヘイム04Fをクリア+LV50 解説 DUEL TERMINAL -シンクロ覚醒!!-(OCG):DUEL TERMINAL -ワームの侵攻!!-(OCG)を収録。 主なカード ブリューナク 収録カードリスト ※レアリティ無表記のカードはノーマル ※OCGの同名パックとの差分は最下段に表記 通常モンスター(??種) 効果モンスター(??種) 儀式モンスター(0?種) 融合モンスター(0?種) シンクロモンスター(0?種) 魔法(??種) 罠(??種) 青眼の白龍(レアリティ未確認) (スーパー) ブラック・マジシャン(レアリティ未確認) (レア) 真紅眼の黒竜(レアリティ未確認) (レア) おジャマ・イエロー(レアリティ未確認) E・HERO ネオス(レアリティ未確認) (レア) バスター・ブレイダー(レアリティ未確認) クリボー(レアリティ未確認) ハネクリボー(レアリティ未確認) サイバー・ドラゴン(レアリティ未確認) (レア) 氷結界の番人ブリズド(レアリティ未確認) (スーパー) ブリザード・ウォリアー(レアリティ未確認) 氷結界の術者(レアリティ未確認) (レア) 霞の谷の雷鳥(レアリティ未確認) 霞の谷の祈祷師(レアリティ未確認) (スーパー) 霞の谷の戦士(レアリティ未確認) フレムベル・ドラグノフ(レアリティ未確認) フレムベル・マジカル(レアリティ未確認) ガード・オブ・フレムベル(レアリティ未確認) X-セイバー アクセル(レアリティ未確認) X-セイバー エアベルン(レアリティ未確認) (スーパー) X-セイバー ウルズ(レアリティ未確認) 総剣司令 ガトムズ(レアリティ未確認) (スーパー) A・O・J クラウソラス(レアリティ未確認) A・O・J ガラドホルグ(レアリティ未確認) A・O・J ルドラ(レアリティ未確認) ワーム・アポカリプス(レアリティ未確認) ワーム・バルサス(レアリティ未確認) ワーム・カルタロス(レアリティ未確認) ワーム・ディミクレス(レアリティ未確認) ワーム・イーロキン(レアリティ未確認) (スーパー) 氷結界の龍 ブリューナク(レアリティ未確認) (ウルトラ) ミスト・ウォーム(レアリティ未確認) (ウルトラ) フレムベル・ウルキサス(レアリティ未確認) (レア) X-セイバー ウルベルム(レアリティ未確認) (ウルトラ) A・O・J カタストル(レアリティ未確認) (レア) クロス・ソウル(レアリティ未確認) 悪魔のくちづけ(レアリティ未確認) 撲滅の使徒(レアリティ未確認) 滅びの爆裂疾風弾(レアリティ未確認) (レア) 黒・魔・導(レアリティ未確認) 黒炎弾(レアリティ未確認) (レア) ラス・オブ・ネオス(レアリティ未確認) (ウルトラ) 機雷化(レアリティ未確認) バーサーカークラッシュ(レアリティ未確認) エヴォリューション・バースト(レアリティ未確認) (レア) 砂塵の大竜巻(レアリティ未確認) 弱体化の仮面(レアリティ未確認) マジック・ジャマー(レアリティ未確認) 援軍(レアリティ未確認) 攻撃の無力化(レアリティ未確認) ニトロ・シンクロン(レアリティ未確認) (レア) デッド・ガードナー(レアリティ未確認) ビッグ・ピース・ゴーレム(レアリティ未確認) (レア) ダーク・リゾネーター(レアリティ未確認) 手錠龍(レアリティ未確認) 疾風の暗黒騎士ガイア(レアリティ未確認) ハーピィ・レディ1(レアリティ未確認) 宝玉獣 サファイア・ペガサス(レアリティ未確認) ジェネクス・コントローラー(レアリティ未確認) (レア) ジェネクス・ウンディーネ(レアリティ未確認) (スーパー) ジェネクス・パワー・プランナー(レアリティ未確認) ジェネクス・サーチャー(レアリティ未確認) ジェネクス・ワーカー(レアリティ未確認) 氷弾使いレイス(レアリティ未確認) 氷結界のロイヤル・ナイト(レアリティ未確認) (スーパー) 氷結界に住む魔酔虫(レアリティ未確認) (レア) 霞の谷の見張り番(レアリティ未確認) ミスト・コンドル(レアリティ未確認) (レア) 霞の谷の風使い(レアリティ未確認) フレムベル・パウン(レアリティ未確認) (レア) フレムベル・アーチャー(レアリティ未確認) フレムベル・デビル(レアリティ未確認) (スーパー) X-セイバー アナペレラ(レアリティ未確認) X-セイバー ガラハド(レアリティ未確認) (レア) X-セイバー パロムロ(レアリティ未確認) X-セイバー パシウル(レアリティ未確認) A・ボム(レアリティ未確認) (レア) A・O・J サウザンド・アームズ(レアリティ未確認) A・O・J アンノウン・クラッシャー(レアリティ未確認) (レア) ワーム・ファルコ(レアリティ未確認) ワーム・グルス(レアリティ未確認) ワーム・ホープ(レアリティ未確認) ワーム・イリダン(レアリティ未確認) (スーパー) ワーム・ジェートリクプス(レアリティ未確認) (レア) ワーム・キング(レアリティ未確認) (ウルトラ) ニトロ・ウォリアー(レアリティ未確認) (スーパー) ゴヨウ・ガーディアン(レアリティ未確認) (スーパー) ハイドロ・ジェネクス(レアリティ未確認) (ウルトラ) X-セイバー ウェイン(レアリティ未確認) (ウルトラ) A・O・J ライト・ゲイザー(レアリティ未確認) (ウルトラ) 痛み分け(レアリティ未確認) フォース(レアリティ未確認) サイクロン(レアリティ未確認) 黒いペンダント(レアリティ未確認) メテオ・ストライク(レアリティ未確認) 道連れ(レアリティ未確認) 誘惑のシャドウ(レアリティ未確認) 粘着落とし穴(レアリティ未確認) 宝玉の祈り(レアリティ未確認) 昇天の黒角笛(レアリティ未確認) ※OCGの同盟パックとの差分 通常モンスター(??種) 効果モンスター(??種) 儀式モンスター(0?種) 融合モンスター(0?種) シンクロモンスター(0?種) 魔法(??種) 罠(??種)
https://w.atwiki.jp/javadsge/pages/4093.html
!DOCTYPE html html head title 統計的方法 /title meta charset="UTF-8" link rel="stylesheet" href="material.min.css" link rel="stylesheet" href="https //storage.googleapis.com/code.getmdl.io/1.0.6/material.indigo-pink.min.css" script src="https //storage.googleapis.com/code.getmdl.io/1.0.6/material.min.js" /script link rel="stylesheet" href="https //fonts.googleapis.com/icon?family=Material+Icons" script src="sample.js" /script script src="jquery.js" /script script src="ols.js" /script /head body onload="sample()" div class="mdl-layout mdl-js-layout mdl-layout--fixed-header" header class="mdl-layout__header" div class="mdl-layout__header-row" span class="mdl-layout-title" 統計的方法 /span /div div class="mdl-layout__tab-bar mdl-js-ripple-effect" a href="#scroll-tab-1" class="mdl-layout__tab is-active" 表 /a a href="#scroll-tab-2" class="mdl-layout__tab" 解答 /a a href="#scroll-tab-3" class="mdl-layout__tab" 出所 /a /div /header div class="mdl-layout__drawer" span class="mdl-layout-title" タイトル /span /div main class="mdl-layout__content" section class="mdl-layout__tab-panel is-active" id="scroll-tab-1" div class="page-content" id="tab2" /div /section section class="mdl-layout__tab-panel" id="scroll-tab-2" div class="page-content" id="memo" /div /section section class="mdl-layout__tab-panel" id="scroll-tab-3" div class="page-content" 13章2節 /div /section /main footer class="mdl-mini-footer" div class="mdl-mini-footer__left-section" div class="mdl-logo" リンク /div ul class="mdl-mini-footer__link-list" li a href="http //osomatsusan.com/" コード /a /li li a href="https //www.google.com/fusiontables/data?docid=1IN8lVhSneZMzzP4bMBG_LDmxZZFrwTik-ukJy0Pr#rows id=1" データ /a /li /ul /div /footer /div /body /html var x1 = new Array(); var x2 = new Array(); var y = new Array(); var x = new Array(); var sol = new Array(); function sample() { var query = "SELECT * FROM " + 1aM2SHx_gjN2Ipbcac0mbqt2Gps47ceaedXMfSUUC ; var encodedQuery = encodeURIComponent(query); var url = [ https //www.googleapis.com/fusiontables/v1/query ]; url.push( ?sql= + encodedQuery); url.push( key=AIzaSyAH0WCAXRIdHdbI-6bnMSN4kVDDZZmGyqY ); url.push( callback=? ); $.ajax({ url url.join( ), dataType jsonp , success function (data) { var rows = data[ rows ]; sx = 0; for (var i in rows) { sx = sx + 1; x1[sx] = rows[i][0] / 1; x2[sx] = rows[i][1] / 1; y[sx] = rows[i][2] / 1; } number = sx; for (s = 1; s number + 1; s++) { x[s] = new Array(); } for (s = 1; s number + 1; s++) { x[s][1] = 1; x[s][2] = x1[s]; x[s][3] = x2[s]; } var sub = new ols(); sub.x = x; sub.y = y; sub.numberx = 3; sub.datanumber = number; sub.mat(); sol = sub.sol; paint(); paintx(); } }); } function paint(){ str1=" table class=\"mdl-data-table mdl-js-data-table mdl-data-table--selectable mdl-shadow--2dp\" "; str2=" th class=\"mdl-data-table__cell--non-numeric\" "; str3=" td class=\"mdl-data-table__cell--non-numeric\" "; str=str1; str=str+" thead "; str=str+" tr "; str=str+" th x1 /th "; str=str+" th x2 /th "; str=str+" th y /th "; str=str+" /tr " str=str+" /thead "; str=str+" tbody "; for(s=1;s number+1;s++){ str=str+" tr "; str=str+" td "+x1[s]+" /td "; str=str+" td "+x2[s]+" /td "; str=str+" td "+y[s]+" /td "; str=str+" /tr "; } str=str+" /table "; $("#tab2").html(str); } function paintx(){ $("#memo").html(sol[1]+","+sol[2]+","+sol[3]); } class ols{ x any=new Array(); y number[]=new Array(); datanumber number; numberx number; sol number[]=new Array(); a any[]=new Array(); b number[]=new Array(); mat() void{ var s number; for(s=1;s this.numberx+1;s++){ this.a[s]=new Array(); } var s1 number; var s2 number; var d1 number; for(s1=1;s1 this.numberx+1;s1++){ for(s2=1;s2 this.numberx+1;s2++){ d1=0; for(s=1;s this.datanumber+1;s++){ d1=d1+this.x[s][s1]*this.x[s][s2]; } this.a[s1][s2]=d1; }} for(s1=1;s1 this.numberx+1;s1++){ d1=0; for(s=1;s this.datanumber+1;s++){ d1=d1+this.x[s][s1]*this.y[s]; } this.b[s1]=d1; } var z number; for(s1=1;s1 this.numberx+1;s1++){ for(s2=1;s2 this.numberx+1;s2++){ z=this.a[s2][s1]/this.a[s1][s1]; if(s1==s2)z=0; for(s=1;s this.numberx+1;s++){ this.a[s2][s]=this.a[s2][s]-z*this.a[s1][s]; } this.b[s2]=this.b[s2]-z*this.b[s1] } } for(s=1;s this.numberx+1;s++){ this.sol[s]=this.b[s]/this.a[s][s]; } } }
https://w.atwiki.jp/javadsge/pages/3026.html
(1)プログラム code (2)グラフ (3)出所 apache (4)メモ 地理 (5)作業記録 11月4日 追加 -
https://w.atwiki.jp/javadsge/pages/3107.html
(1)プログラム code (2)グラフ (3)出所 apache (4)メモ 1.5 重回帰 (5)作業記録 11月4日 追加 -
https://w.atwiki.jp/javadsge/pages/1293.html
!DOCTYPE html html lang="ja" head meta charset="utf-8" / script src="sample.js" /script script src="ols.js" /script script src="jquery.js" /script script src="svgpaint.js" /script title 回帰 /title /head body onload="sample()" p id="memo" /p /body /html var x = new Array(); var y = new Array(); var sol = new Array(); var datanumber,numberx; function sample(){ var data; data = new XMLHttpRequest(); data.open("GET", "data.txt", false); data.send(null); var str=data.responseText; var g=str.split("\n"); datanumber=g.length; for (s = 1; s datanumber+1; s++) { x[s] = new Array(); } for(s=0;s datanumber;s++){ var g1=g[s].split(","); y[s+1]=g1[0]; x[s + 1][1] = 1; x[s+1][2]=g1[1]; } numberx = 2; var sub = new ols(); sub.x = x; sub.y = y; sub.numberx =numberx; sub.datanumber =datanumber; sub.mat(); sol= sub.sol; paint(); } function paint(){ var sub6=new svgpaint(); str=" svg width=\"700\" height=\"700\" "; for(s=1;s datanumber+1;s++){ z1 = sol[1] + sol[2] * x[s][2]; str=str+sub6.circle(x[s][2],500-y[s],5,"blue"); str=str+sub6.circle(x[s][2],500-z1,5,"red"); } str=str+sub6.line(0,0,0,500,2,"black"); str=str+sub6.line(0,500,500,500,2,"black"); str=str+" /svg "; $("#memo").html(str); } class ols{ x any=new Array(); y number[]=new Array(); datanumber number; numberx number; sol number[]=new Array(); a any[]=new Array(); b number[]=new Array(); mat() void{ var s number; for(s=1;s this.numberx+1;s++){ this.a[s]=new Array(); } var s1 number; var s2 number; var d1 number; for(s1=1;s1 this.numberx+1;s1++){ for(s2=1;s2 this.numberx+1;s2++){ d1=0; for(s=1;s this.datanumber+1;s++){ d1=d1+this.x[s][s1]*this.x[s][s2]; } this.a[s1][s2]=d1; }} for(s1=1;s1 this.numberx+1;s1++){ d1=0; for(s=1;s this.datanumber+1;s++){ d1=d1+this.x[s][s1]*this.y[s]; } this.b[s1]=d1; } var z number; for(s1=1;s1 this.numberx+1;s1++){ for(s2=1;s2 this.numberx+1;s2++){ z=this.a[s2][s1]/this.a[s1][s1]; if(s1==s2)z=0; for(s=1;s this.numberx+1;s++){ this.a[s2][s]=this.a[s2][s]-z*this.a[s1][s]; } this.b[s2]=this.b[s2]-z*this.b[s1] } } for(s=1;s this.numberx+1;s++){ this.sol[s]=this.b[s]/this.a[s][s]; } } }
https://w.atwiki.jp/ce00582/pages/1103.html
import java.awt.*; class prob{ double[] ex(){ double p; double mu; double sig; double yy; double th[]= new double[101]; int s; mu=0; sig=0.39; for (s=1;s 101;s++){ p=0.01*s-0.005; yy=seeky(p,mu,sig); th[s]=Math.exp(yy); } return th; } static double f(double x,double mu,double sig){ double pi,x1,x2,x3,fx; pi = 3.1415; x1=-Math.pow(x - mu,2) / (2*Math.pow(sig,2)); x2=Math.exp(x1); x3=sig*Math.pow(2*pi,0.5); fx=x2/x3; return fx; } static double g(double y,double mu, double sig){ double gx,h,x; int n,t; gx=0; h=0.001; t=(int)(y/h); for (n=-2000;n t;n++){ x=n*h; gx=gx+f(x,mu,sig)*h; } return gx; } static double seeky(double p,double mu,double sig){ double g1,g2,y1,y2,y3; int t; y1=0.4; y2=-0.2; g1=g(y1,mu,sig); t=0; while(t 100){ g2=g(y2,mu,sig); y3=y2+(p-g2)*(y2-y1)/(g2-g1); y1=y2; y2=y3; g1=g2; if (Math.pow(p-g2,2) 0.0001)t=1000; t=t+1; } return y2; } } class Can extends Frame { public static void main(String[] args) { new CanvasTest(); } Can() { super("CanvasTest"); setSize(200, 100); setLayout(new BorderLayout()); MyCanvas mc1 = new MyCanvas(); add(mc1, BorderLayout.CENTER); show(); } } class MyCanvas extends Canvas { public void paint(Graphics g) { int thp[]=new int[101]; double th[]=new double[101]; int s; prob ton =new prob(); th=ton.ex(); for (s=1;s 101;s++){ thp[s]=(int)(30*th[s]); } for (s=1;s 90;s++){ g.drawLine(s,thp[s], s+1,thp[s+1]); } } }